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ABSTRACT 

A previous analysis of the Reuven Ramaty High Energy Solar Spectroscopic 
Imager (RHESSI) observation of GRB 021206 found that the gamma-ray flux 
was 80±20% polarized. We re-examine this data and find no signal that can be 
interpreted as due to polarization. First, we find that the number of scattering 
events suitable for measuring polarization - having been scattered from one 

^ ■ detector to another, with a count produced in both - is considerably lower 

than estimated by CB03, by a factor of 10 (830 ± 150, vs. 9840 ± 96). The 
> ■ signal-to-noise of the data-set is thus too low to produce a detection, even from a 

100% polarized source. Nonetheless, we develop a polarization-detection analysis 
limited in sensitivity only by Poisson noise, which does not require a space-craft 

^ ' mass model to detect polarization, as in CB03. We find no signal which might 

be interpreted as due to polarization of GRB 021206. Separately, we reproduce 

■ the CB03 signal and show that it is not due to polarization. Rather, the CB03 
signal is consistent with the previously-neglected systematic uncertainty in the 

q "null lightcurve" used for detection. Due to the low signal-to-noise ratio of 

the RHESSI data, our Poisson noise-limited analysis results in an upper limit 
consistent with 100%-polarization of the gamma-ray flux from GRB 021206. 

■ Thus, no observational constraint on the polarization of GRB 021206 can be 
^ ■ derived from these data. 

.5?: 

1. Introduction 

Gamma-ray bursts (GRBs) are associated with the deaths of massive stars, as seen 
first in their association with blue galaxies (Bloom et al. 1999; Bloom et al. 2002a); 
the appearance of supernovae explosion (SNe) -like lightcurves at late times (Bloom 
et al. 2002b); and the recent convincing observation of a SNe spectrum associated with a 
GRB afterglow (Hjorth et al. 2003; Stanek et al. 2003). However, what differentiates GRBs 
from the 10 4 more frequently observed core-collapse, or type-II, SNe remains uncertain. 
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Recently, it was reported (Coburn & Boggs 2003, CB03 hereafter) that 80±20% 
polarization in the prompt, or burst, gamma-rays of GRB021206 was detected using the 
Spectroscopic Imager (Smith et al. 2002) on board the RHESSI satellite (Lin et al. 2002). 
Although not designed for gamma-ray polarization measurements of GRBs, the multiple 
detectors of RHESSI were used to search for simultaneous events in two detectors due to 
Compton scattering, which then give a preferred scattering direction projected on the sky. 
An excess of scattered events in a particular direction on the sky was interpreted as due to 
angle-dependent scattering associated with intrinsic polarization of the gamma-ray photons. 

Such a high polarization fraction in the gamma-rays may require a large-scale persistent 
magnetic field - of a strength observed thus far only from neutron stars, in particular 
the magnetars which are believed to account for <10% of the observed NS population. 
If this property of GRBs is confirmed in future observations, then polarization would 
strongly constrain GRB production models, and may provide important input physics 
to understanding the generalized SNe phenomena. The reported detection has led to 
the wide discussion of mechanisms for producing the high polarization, which would 
constrain the emission and progenitor models (Lyutikov et al. 2003; Granot 2003; Eichler 
& Levinson 2003; Nakar et al. 2003; Lazzati et al. 2003; Matsumiya & Ioka 2003), and has 
prompted discussion of detector development to better observe gamma-ray polarization 
(Bloser et al. 2003). 

The analysis of CB03 used a mass-model for the RHESSI satellite to calculate the 
satellite response to an unpolarized beam, in order to calculate the unpolarized lightcurve. 
It is first needed to calculate the "null" lightcurve (that is, the lightcurve of an unpolarized 
GRB) to demonstrate detection. It is also needed to correct the detected signal, and 
dominates the error bar in the measurement of 80±20% polarization. 

We describe the observation in §1.1, and give an overview of the analysis by CB03 
in §1.2. To produce the polarization measurement, it is necessary that photons scattered 
between two detectors be observed in the data, and in §2 we estimate the number of 
observed scattered double-count events, finding it to be a factor of 10 lower than the 
previous estimate. In §3, we examine the distribution of single counts and double counts 
as a function of the instantaneous position of the detectors where the counts are detected, 
and find that, while both single- and double- counts exhibit strong dependencies on this 
position, their ratio does not, suggesting the absence of a signal due to polarization. We 
nonetheless describe in §4 a polarization detection analysis suitable for RHESSI data which 
is limited in its sensitivity by Poisson statistics (that is, not dependent upon the space-craft 
mass model) and apply it to the observation of GRB 021206, deriving a correction which 
accounts for the angle- dependent Klein-Nishina cross section in §4.2. The results of the 
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polarization analysis are discussed in §4.3, finding no evidence for a signal which could be 
due to polarized gamma-ray photons from GRB 021206. In §5, we reproduce the modulation 
reported by CB03, finding that the reported signal is not observed using the new method, 
which does not rely upon the Monte Carlo radiative transfer through the RHESSI mass 
model. In §8 we summarize our results, compare them with those of CB03 and conclude. 

1.1. Observation of GRB 021206 with RHESSI 

The observation of GRB 021206 is described by CB03. We recount pertinent 
information. The RHESSI satellite is a gamma-ray imager primarily for solar observations 
using a rotation modulated collimator. It is in a low-Earth orbit, with its imaging (z) 
axis staring within a few arc-minutes of the solar center. The entire satellite rotates 
approximately around this z axis clockwise from north, with a period of ~4 seconds. 

The RHESSI spectrometer detectors are an array of 9 separate, effectively identical 
detectors, each a cylinder roughly 7.1 cm diameter. Each detector is electronically 
separated into two segments ("front" and "rear") which have slightly different photon 
energy responses. A diagram of the detectors and their orientation within the spacecraft 
detector plane is given in (McConnell et al. 2002). We give a similar diagram in Fig. 1, and 
the mass-model coordinate centers of the detectors in Table 1. 

The approximate start time of GRB021206 measured with RHESSI was 2002-Dec-06 
22:49:16 UT. The best gamma-ray localization of GRB021206 was found (Hurley et al. 2003) 
to be an error ellipse centered at R.A.=240.195 deg (16h00m46.8s), dec. =-9. 710 deg 
(— 09d42m36s)(J2000), with a major axis of 20.4'and a minor axis of 0.53', position angle of 
-18 deg (relative to the positive direction of declination). We adopt the position of a related 
radio transient (Frail 2003). 

At the time of the transient, the solar location was about R.A.=16h53m20.15s, 
dec.=-22d32m47.7s. The GRB was localized at an angle #grb = 45 deg clockwise from N 
of the sun, and offset from the direction of solar center by 17.98 deg. Counts detected with 
the RHESSI spectrometer are timestamped with resolution of 2 -20 s, or 1 binary micro-sec 
(b/xs); the detector segment; and photon energy. 

1.2. Approach of CB03 

We refer the reader to CB03 for details of the analysis used to detect the polarization 
signal. We describe the pertinent approach here. 
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A double-count scattering event takes place when a gamma-ray photon makes a single 
scatter off of detector i, producing a count in that detector, then leaving that detector to 
be partially or fully absorbed in a second detector j, producing a simultaneous count in 
detector j. A double-count coincidence event, which appears identical to a double-count 
scattering event, takes place when two unrelated photons produce simultaneous counts in 
two different detectors. A single-count event takes place when a single photon is fully or 
partially absorbed in a detector i, but is not followed by a detection in another detector 
within some period of time AT as measured by the on-board clock. 

When a double-count event is registered, the position angle on the sky (6) of the line 
joining the centers of the detector pair is determined. From this, a lightcurve of double 
count events as a function of 9 was produced. 

CB03 produced a null double-count event lightcurve to be subtracted from the observed 
lightcurve by Monte Carlo simulation (S. Boggs, W. Coburn, priv. comm.). They simulated 
~18,000,000 photons propagating through the spacecraft from the direction of the GRB, to 
produce a library of double-event scatters between detectors. 

Following this subtraction, a sinusoidal modulation in the difference lightcurve was 
found, which was significantly different from a constant value. Based on this sinusoidal 
residual, CB03 concluded that the incoming gamma-ray beam was polarized. Using the 
scattering fractions found with the GEANT mass model, the magnitude of the modulation 
was corrected to find the intrinsic polarization magnitude of 80±20%. The uncertainty is 
dominated by uncertainty in the GEANT mass model. 

2. Proportion of Scattering Events and Backgrounds 

To estimate the sensitivity of the RHESSI spectrometer to gamma-ray polarization, 
it is critical to establish the presence and estimate the number of two-detector scattering 
events in the data, and to minimize, and also estimate, the level of contaminating 
two-detector backgrounds. 

It is possible to analyze the RHESSI data for evidence of a polarization signature 
without investigating these effects. In most but not all cases (see §6), double-count 
event backgrounds will dilute the angular modulation signal but will not generate an 
angularly-modulated signal of their own, and thus will not mimic the effects of genuine 
gamma-ray polarization. However, even when modulation can be detected in this fashion, 
estimating the number of two-detector scatters and associated backgrounds still allows a 
valuable consistency check of the analysis to be made. Moreover, converting a detection of 
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modulation to an estimate of the polarization signal strength requires an accurate estimate 
of the signal and background. 

CB03 present estimates of the number of two-detector scatter events and backgrounds 
in their "Methods" section. They find 9840 ± 96 two-detector scatters, and estimate a 
background of 4488 ± 72 two-detector coincidences and 588 ± 24 "background scatter" 
events, for a total of 14916 double-count events detected. The meaning of the "background 
scatter" category is not defined but would seem to refer to scattering events attributable 
to background sources other than GRB 021206 itself. They also state that the number 
of two-detector scatter events is "roughly 10% of the total 0.15-2.0 MeV light-curve 
events," and that this proportion is in agreement with the results from their Monte Carlo 
simulations; the nature and degree of this agreement is not quantified. 

In this section we make a model-independent investigation of the fraction of two- 
detector scatter events in the RHESSI data for GRB 021206. In particular, we determine 
the number of two-detector events which are due to events of each of the following types: 
(a) scattered events, which produce two counts in two different detectors from a single 
scattered photon; (b) the coincident arrival of two unrelated photons in two different 
detectors; (c) scattered events, which produce two events in the same detector; and (d) 
events which take place in different detectors due to a detector effect. For the polarization 
measurement, events of type (a) permit measurement of polarization; type (b) is the 
irreducible background; and types (c) and (d) are background which, in principle, can be 
identified by their non-Poissonian nature and removed from the data. 

2.1. Relative Timing Accuracy of RHESSI Detectors 

As a first step, we examine the timing properties of RHESSI spectrometer data to 
establish an appropriate definition for selection of "simultaneous" events. We find that it is 
necessary to accept events within a time window of 5 b/is to avoid discarding simultaneous 
events. 

Fig. 2a is a histogram of wait times between timestamps of consecutive counts (AT); for 
a pure Poisson process - not the present case, because of the strong variations in the GRB 
mean count rate - this histogram will exhibit the classical exponential form. We observe a 
clear non-Poissonian excess near AT pa 2 b/is, that is, at AT greater than zero. Since the 
light-travel time for 2 b/is is 570 m, far larger than the separations between detectors, this 
cannot be a physical effect, but rather must be due to behavior of the detector electronics, 
with simultaneous events not receiving identical timestamps. 
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In Fig. 2b we have removed the 6732 event pairs which take place either simultaneously 
or consecutively (with no intervening counts) in the same detector; the front and rear 
segments of a detector are considered to be the same detector for these purposes. Note 
that double-count, same-detector events are not useful for polarization measurement. 
While these events are a significant fraction of those with AT < 5, we still observe a 
non-Poissonian excess of counts at AT=l-4b/is; in addition, we can see non-Poissonian 
excesses at AT=8, 12, 16, 20, 24, 28 and 32 hfis 2 . These excesses indicate that the 
spacecraft timestamps are subject to systematic uncertainties at the AT ~ 4b/is level, 
and possibly even at the 32b/is level. Thus, we treat event pairs with AT < 5b/zs as 
simultaneous in our subsequent analyses. 

As a consequence of this definition, if we observe more than two consecutive counts 
within 5b/is of each other, their simultaneity makes it impossible to extract a unique 
scattering event from the group. We therefore exclude all events within such multi-event 
groupings from our analysis. 



To evaluate the degree to which the counts statistics are Poissonian, we bin the data 
into a lightcurve with time resolution AT=5 b/is. 

To account for the highly variable nature of the GRB light curve, we estimate the 
countrate in St = 0.005 s intervals, treating the count rate over each interval % as constant 
for purposes of estimating the Poissonian mean rate /Xj. For the entire selected portion of 
the GRB light curve, we find the number Aj j0 b s of 5 b/xs bins containing j counts to be 
AW =969787, A lj0bs =72510, A 2>obs =5666, and A >2i0bs =481 bins, respectively. 

Estimating the Poisson mean rate /ij = A ljObSi j/A 0iObSi j in each interval from the number 
of single-event and zero-event time bins, we calculate the number of 5 b/is time bins within 
that interval that will contain 2 (A 2 ) or >2 (A >2 ) events, respectively, due to Poisson 
fluctuations: 



2.2. 



Double-Count Events Due to Coincidence 



= ^E^ 
AT^ 2 



(1) 



N- 



>2 



St 
AT 



El-(l + /^ + y)e 



(2) 



2 The excesses are even more apparent if one looks only at events detected in the rear segments. The 
origin of these excess counts at 4byus is related to a detector deadtime effect (D. Smith, priv. comm.) 
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Summing over all St intervals, we find N 2 = 3760 ± 30 bins. The remaining excess of 
(5666 — 3760) = 1906 ± 81 double-event bins are not due to coincidence (event type b), and 
must be due to non-Poisson processes (event types a, c, or d). 

A disproportionate fraction of these arise from double-count events within a single 
detector (event type c). Of the iV 2)0 bs — 5666 events, 1313 are due to simultaneous events 
in the same detector, compared to the 1/9 x 5666 630 events expected from coincidence. 

As a result, when we require the two events in each double-event bin to rise from 
different detectors, we find an observed number of A 2i0bs ,dd = 4983 ± 83 double-event, 
different-detector bins compared to an estimated number of A 2jdd = 8/9 x N 2 = 3342 ± 27 
due to Poisson coincidence alone. This implies that 1641 ±87 double-event, different-detector 
bins may be due to scattering or unidentified detector effects (event types a or d). 

2.3. Double- and Multi-Count Events Due to an Unknown Process 

While we find the number of bins with > 2 events due to coincidence is A >2 =159±2 
bins, we observe A >2i0 bs.=481 bins, for an excess of 322±22 bins, which must be due to 
non-Poisson processes (event types a, c or d). 

If these non-Poisson events were due exclusively to scattering processes, then the 
number of double-count events will be a fraction / of single events A 2iScat terin g = f Ni, the 
number of triple scattering events A 3iScattcring = f 2 N ± , and number of j- scattering events 
Aj iSca ttering = P' 1 N ± . Thus, the ratio r = A 2iScattcring /A >2iScattcring = l./(£ i=1 /'). Taking 
/ = A 2i0bs ./Ai i0bs . = 1906/85387=(2.2±0.1)xl0~ 2 , we find r = 44 ± 2, while we observe 
^obs. = A 2j0 bs./A >2i0 bs. = 1906/322 = 5.9±0.5 (here, we have not corrected for double-count 
events due to scatters within a detector) 3 . 

Thus, there are a factor 7.8±0.8 more A >2 o b s . events relative to the number of A 2 
events than can be explained by scattering events. We compared the distribution of 
detectors and photon energy for counts included among the observed A >2 events, and found 
they do not differ significantly from the same distributions of all the GRB data. Specifically, 
the multi-count events typically occur in more than one detector. 

We performed an identical analysis on data from during a background period 



3 The cross section for Compton scattering increases by a factor of ~2 between 1 MeV and 150 keV. In 
a worse-case, the ^scattering = (2/)/iVi and iV N>2 , S cattcring = {2f) N ~ 1 fN 1 , and / = 22.7 - still discrepant 
with the observations. 
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(43180 counts in 24 s), and found non-Poissonian N >2 bins present in those data as 
well. The background period should have produced iV >2ibg = 0.2 bins (on average), but 
A^2,b g ,obs.— 1013. Also during this period, N 2j b g = 105 ± 1, while iV 2i b gi obs.=3827 bins were 
observed. Therefore, these non-Poissonian excess events are not associated with the GRB. 

We attribute this non-Poissonian excess to an unknown process - perhaps particle 
background, or noise in the detector electronics or instrumentation - which we are unable 
to precisely model. Nonetheless, based on the background period, the non-Poisson process 
produces a ratio of counts r b s .,b g = (3827 — 105)/(1013 — 0.2) = 3.7 ± 0.1, significantly 
below the value of 5.9±0.5 observed during the burst. 

It seems reasonable that the same unknown process can also produce double-count 
events, contributing to the double-count event background; however, the process is clearly 
non-Poissonian, and we are unable to produce a demonstrably reliable estimate of the 
number of such events 4 . Nonetheless, we note the effect and its uncertain contribution to 
the irreducible background as an event type d. 

2.4. Summary of Corrections for Double-Count Event Rate 

From analysis in this section, we conclude: 

1. counts with time-tags different by < 5 b/zs should be regarded as simultaneous; 

2. double-count events within a single detector contribute significantly, adding to the 
background of double-count events, but are easily identified in the data stream, and 
should be replaced with a single count at the earlier time-stamp; 

3. An unknown effect produces multi-count groups (> 2) of simultaneous counts in 
different detectors in excess of the number expected from coincidence, and which 
cannot be attributed to a scattering process, both during the GRB and during an 
earlier background period. These groups can be identified in the data-stream as 
occurring in a period <5b/is; since they cannot be used for measuring polarization 
- as which two of the > 2 events were the sequentially first cannot be determined - 
these events should be removed prior to analysis; 



4 If, however, we assume that r bs. is constant between the background and the GRB observation, then the 
number of double-count events during the GRB which could be due to type a events is 1906±81 - 3.7±0.1x 
(322±22) = 715±120 counts. 
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4. The same non-Poissonian, non-scattering effect likely also produces double-count 
events; we can estimate the number of double-count events based on the ratio 
r obs. = (-^2,obs. — N 2 ) / (N >2j0 bs. — ^>2) (correcting each value for the coincidence 
rate) observed during the background observation. However since we cannot model 
the non-Poissonian process which is responsible for this ratio, the ratio has an 
unquantifiable systematic uncertainty 

Based on this analysis, we sequentially apply the following selections to the GRB data 
to obtain double-count events using un-binned data: 

• Replace all double-count events in the same detector with single-counts, at the earlier 
time stamp (83300 counts remaining afterwards). This removes "echo" counts in the 
detector (Smith 1998). 

• Remove groups of >2 events which take place in <5 b/is, (iV >2i0bs . =719 such 
multi-count events, which is in excess of the iV >2 =474 such events expected from 
Poisson statistics, leaving 81034 counts). 

Following these selections, we have A^ 2 , o bs.=8230 double-count events remaining. This 
is fewer than the iV 2)0bs . = 14916 found by CB03. We list the contributions of "signal" and 
the irreducible background to these totals in Table 2. 

To estimate the contribution to the irreducible background caused by coincidence of 
two unrelated photons in two different detectors, we use a similar approach as in § 2.2. 
We use an integration time of St =0.005 sec, and find the number of single counts which 
have another which follows within AT = 5b/is. To do so, we calculate N n , the number of 
unbinned events which are a collection of n counts within a period AT: 

Nn = Y. (fri)! e M-mm - (n - 1)) 

where JVj is the number of detected counts in bin i, and fa = NiAT/St. Using this we find 
iV 2 =6640±80, after correcting by a factor 8/9 to account for the requirement that the counts 
be in two different detectors. This is greater than that found by CB03 (iV 2 =4488±72). If 
we were to assume a constant countrate throughout the burst we obtain N 2 = 5380 - lower 
than what we find for a variable intensity burst, but above that of CB03. 

To estimate the contribution to the irreducible background due to the 
"unknown" process, we performed the same analysis as in §2.3, but to the 
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unbinned data, obtaining r obs =3.1±0.1 from the background period, giving us 
A^bkg = r ohs , x (N >2>obs , - N >2 ) = 3.1 ± 0.1(719 - 474) = 760±110 double-count 
events due to the unknown process (this value has an unquantifiable systematic 
uncertainty). 

The result of this is that we find significantly fewer non-coincidence double-count 
events than found by CB03 (830±150 double-count events, vs. 9840±96 by CB03). 

We use these single and double-count event lists for analyses in §2.5, §3 & §4. 
2.5. A Simpler Method of Determining the Double-Count Event Rate due to 



There is a method to estimate the fraction of double-count events which are due to 
scattering between detectors considerably more straightforward than attempting to identify 
all the backgrounds, as we do above. 

First, we note that for each detector the total solid angle subtended by secondary 
detectors varies (see Fig. 1). In particular, the relative solid angle seen by the inner 
detectors (numbers 1, 2, and 7) is significantly larger than that seen by the outer detectors 
(numbers 3-6, 8, and 9), and if scattering events represent a significant fraction of the total 
then this should be reflected in the number of double-count events involving each detector. 

For the analysis in this sub-section only, we exclude detector 2, due to its different 
electrical set-up. We calculate, for each detector, the total physical solid angle (f^) 
subtended by the other seven detectors, except where the lines of sight to the other 
detectors are blocked by detector 2. We also examined the total number of single-counts 
in each detector Jj observed during the GRB - if all detectors were equally sensitive 
(including shadowing by the spacecraft), they should have equal values of /«, but they do 
not, indicating different sensitivities during the GRB. We use Ii as a measure of the detector 
responses (see Table 1 for Qi and If, note that we use the counts selected as described in 



We then determine iV 2i0 bs(^) ~~ the number of double-count events in which one of the 
counts appears in detector i. We model this observed event rate with the function: 



where A^tot is held fixed at the observed number of counts for the double-events in detectors 



Scattering Between Detectors 



§2-4). 




(3) 
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1 and 3-9 (10906). Here, the coefficient on / is due to the scattered events involving two 
detectors. The coefficient on 1 — / is the random background. The value / can be measured 
because the values of Qj are different for different detectors. Our analysis ignores two 
effects: the angular dependence of the Klein-Nishina cross section for scattering, and the 
absorption/scattering effects of passive material between detector units, which will decrease 
the effective Jlj by varying amounts, depending on the distance between detectors. 

Using a x 2 minimization fitting technique, we obtain an acceptable fit to the data 
(x 2 =9.28 for 7 degrees of freedom - dof), for / = 0.11±0.03. The best fit and residuals, 
including the proportion of double events associated with scattering and coincidence, 
respectively, for each detector, are shown in Figure 4. Therefore, we find 11±3% of the 
double-count events are due to scattering between detectors. 

Comparing this with the value we found by attempting to eliminate all the reducible 
background sources we find (Il±3%x8230±90 double-count events) =910±250 double-count 
events due to scattering between detectors (where here use double-count events in 
all detectors, including detector 2). This is in agreement with the value of 830±130 
double-count events we found from the detailed analysis of the different background 
contributions. 



3. Distribution of Observed Counts with Detector Angle 

Fig. 5a presents the number of single counts S as a function of detector angle 9, where 
for each event 9 indicates the position of the detector at the time of the event, relative to 
the center of the detector array in a non-rotating celestial frame (9 = is celestial north). 

We see that there is a clear bias for single events to occur towards one side of the 
spacecraft. This could be due to a number of effects. For example, different detector 
sensitivities produce uneven exposure as a function of time over the range of 9 which, 
combined with the GRB intensity variability (Fig. 3), can produce such a bias. Also, 
Earth-scattered gamma-rays - which may constitute a significant fraction of the total events 
- will preferentially illuminate the side of the spacecraft towards the Earth (see McConnell 
et al. 2002 for discussion). 

The increased rate of single events towards one side of the spacecraft will produce 
an increased rate of double events on the same side through the increased coincidence 
rate (it will also produce an increased rate of double events through scattering). Fig. 5b 
presents the distribution D{9) for the double events, where 2 values of 9 are used (one for 
each detector involved in the double). D{9) shows the same bias as the S{9) distribution. 
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Indeed, the 9 angular variations in the single and double count rates are highly correlated. 
This can be seen in Fig. 5c, which shows that although the xt values for a constant fit to 
the variations of S(9) and D(9) separately are both very high (giving Prob. < 10 -50 in both 
cases), the xl value for their ratio, R{9) = D{9) / S{9), is close to unity, and acceptable for 
a constant fit (Prob. = 0.54). 

If the double count events were due entirely to coincidences of independent single-count 
events we would expect that D{9) oc S(6). However, if a significant fraction of double-count 
events are due to scattering, and polarization of gamma-ray photons produces scatters 
toward a preferred direction in 9 and 9 + ir, as suggested by the results of CB03, then 
we should observe a D{9) that was significantly different than S{9). Our result here 
conflicts with this, suggesting instead that the variability observed by CB03 is not due 
to polarization-effected scattering, but due to systematic uncertainty in the calculation of 
their "null" lightcurve by MC simulation. 

As a result, the coincidence rate for any given pair of detectors will increase when 
they are both on the preferred side of the spacecraft, and this may allow the derived 
"scatter angles" for the coincidence events to pick out a preferred direction on the sky. 
In the next section we describe and perform a polarization analysis which is insensitive 
to contamination from angular structure in the coincidence events, which is clearly the 
dominant effect in the observed double-count event variability. 

4. A Mass-Model Independent Polarization Analysis 

We describe in this and the following sections a method for detecting polarization 
which does not rely upon the mass model of RHESSI . 

The observed countrate in a detector i is: 

Ni(t) = f(t)Mt) 

where f(t) is the time dependent flux (ergcm' 2 s _1 ) and Ai(t) is the detector response 
(counts erg -1 ), which we allow to be time-dependent due to the rotation of RHESSI with 
respect to the sky. We assume that the probability that a count: (1) is a single-scattering 
event in detector i; (2) escapes detector i; and (3) is then partially or fully absorbed in a 
second detector j, is a time-independent factor B it j, dependent only on the materials and 
geometry of the RHESSI spectrometer. 

The number of double-count events due to unpolarized photons between detectors 
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i and j at a time t when the centers of those two detectors produce a line which has a 
position angle 9 i: j on the sky (relative to the north celestial pole) can be found: 



.Y,, ; (/) .Y,( /)/;,,, • .Y,-( /)/>;,• (4) 

If the beam is polarized, this produces a polarization pattern I p (9ij) where &i,j(t) = wt + faj 
is now the angle the detectors make on the sky as a function of time, and u is the rotational 
frequency of the RHESSI detector plane, and faj is the angle the detector pair makes 
on the sky at time t = 0. 5 



Ni,MM) = Ntfuit)) l p (0ij) Bij + Nj(eij(t)) I !,[(>,:, ) B jti (5) 

where we use the polarization intensity pattern: 

1 + p cos(2(fl - 9 P )) 
Ipifi) = —p (6) 

(This differs by ir/2 from a different, oft-used definition). Now, I(O it j) = I(ir + 9 i: j) = I(0j :i ); 
and Bij = Bj :i for the identical detector geometries which are used by RHESSI. Thus: 

W = me^t)) + Njie^mBijiMj) (7) 

Each side of this equation can be summed over all detector combinations, keeping data 
binned as a function of 9; then 

e<9i j <e+A g e<e i3 <e+ a 

£ .V,,("„) E [N l {9 h3 {t)) + N 3 {d t ,{t))]B l , (8) 

i,j>i i,j>i 

which we can then re-write as a ratio, which is a function of sky-angle only: 



m = o< e <iik — — (9) 



5 We do not correct the polarization angle measured at the detector, which is projected relative to the 
North Celestial Pole, to the angle relative to the north from the GRB position in the sky. The difference 
varies as a function of polarization angle by up to 3 deg systematically. This is <20% of our 15 deg bins, and 
will not affect or conclusions. 
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l+pcos(2(g-g p )) 

1+p [ J 

If p = (unpolarized photons) then i?(#) is a constant, with no dependence on 9. 
However, if p ^ 0, then the observed function R{9) will have an angular dependence, as 
I p (9). 

The only contribution to the uncertainty in R(9) is photon counting statistics. More 
specifically, R{9) is independent of the mass model of the spacecraft, and the uncertainty 
of detection of polarization using R(9) does not depend on systematic uncertainties in the 
space-craft mass model as it does in CB03. However, it should be noted that correcting a 
detected signal to find the intrinsic polarization magnitude (that is, turning the detection 
into a quantitative measurement) would still require using the mass-model. 

In practice, the function R(9) is constructed from a discrete sum of double-count events 
and a discrete sum of single-event counts, as we describe in the following subsection. 

Here, we note an additional correction. As discussed in §2.5, the total number of 
single events observed during the GRB Jj varies with detector i, which implies different 
sensitivities to events. If we take the relative sensitivity /; = h/{h) to be independent of 
the detector geometry (which are assumed identical), then Eq. 4 becomes: 



N i j(t)=N i (t)f i B iJ + N i (t)f i B iti (11) 

Here and above, B iy j is a relative solid angle subtended by detector j from detector i. This 
correction affects Eq. 9, which becomes: 



,0<0i, '/ • A„ 



m = ^<,,,ix; , j( t (12) 



4.1. Data Analysis for GRB 021206 

We downloaded data from the public RHESSI archive 6 , which also makes available 
data analysis software we used - the Solar Soft Ware (SSW) system written in IDL. We 
extract counts from photons with energy between 0.15 and 2 MeV during the observation 



6 http://hesperia. gsfc.nasa.gov 
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of the GRB roughly corresponding in time to those analysed by CB03 (see Fig. 3). We find 
85387 such counts in the 5.0 sec interval we analyze, chosen to coincide with that analysed 
by CB03. 

The background number of counts, measured during an earlier observational period, is 
8636±42 counts during the 5.0 second period. 

We find 8230 double-count events (§2.4). To perform the analysis, we use a discrete 
summation for the lightcurve. To construct the 36 independent double-event light-curves 

Nij(0ij), we take: 

0i,j(t k ) = (ut k + fa) mod (tt) (13) 

where t k is the time of the earliest time-stamp of each double-count event. The values of 
(f>ij are relative angles of the line joining the centers of detectors % and j, and the centers of 
detectors 1 and 2, except for which we set equal to the angle on the sky at time t — 
clockwise of N made by the center-line joining detectors 1 and 2, and oo = 27r/(4.09s) (the 
rotational frequency of the spacecraft at the time of the observation). 

The angles 0jj were found using the coordinate positions for each detector as listed in 
Table 1. We divide the 180 degree range into 12 equal non-overlapping segments, beginning 
with 6 n = n x 15deg, with n — [0, 11]. We produce the summed double-event lightcurve 
dependent upon sky angle D(6 n ): 

£(0„) = £Wo-(Mt*)) 
k 

Where N it j takes the value of 1, for a double-event involving detectors % and j with 
time-stamp t k . The angles are binned as 6 n < 9 it j < 6 n + 6$. 

Thus, for each double-event at time t k involving detectors % and j, we calculate the 
angle on the sky from Eq. 13, and add a count to the appropriate 6 bin. This is shown in 
Fig. 5e. 

We then construct the denominator from the single-event light curve: 

S(*n)=E E fjNiMM) 
k ifc,j!=ifc 

where N ih takes the value of 1, fj is described in §4, and for every event k in detector i k , 



-16- 



we calculate the angle on the sky for every detector pair which includes i k 7 . Thus, for each 
count k, detected at time t k in detector i k we find the eight angles on the sky 9 ik j produced 
by the % k and jth detector, and add one count in each of the eight bins, including angle 
bins in which the same count falls more than once due to co-aligned detector pairs (such as 
detector pair 4 & 7 and pair 4 & 3). This is shown in Fig. 5d. 

We produce R(9 n ) = D(9 n )/S(9 n ), which is shown in Fig. 5f. We discuss this Figure 
in detail in §4.3, after first discussing the effect of the Klein-Nishina cross section on the 
detection in §4.2. 

The uncertainty in R(9 n ) is set by the photon counting statistics of D and S. As a 
minor correction, it is the number of unique counts in each bin of S which contributes to 
the fractional uncertainty of S. Thus, if there are 40000 counts in a particular angle bin of 
S, but only 30000 of those are unique, then the contribution to the fractional uncertainty 
of R{9) is 1./V30000. Even so, the uncertainty is dominated by the smaller number of 
double-event counts, which are in the range 610-770 in each bin. 



The GRB was offset from the pointing axis, which is perpendicular to the plane of 
detectors, by 17.98 deg. Assuming a thin detector plane 8 , photons which first scatter in a 
detector on the side of the detector plane toward the GRB can scatter through an angle 
(to the second detector) of up to #kn=90 deg+18 deg=108 deg; while photons which first 
scatter on the side of the detector plane away from the GRB could scatter through an angle 
of up to 9 KN =90-18=72 deg. 

The difference in scattering geometry produces a varying scattering cross section for 
gamma-rays in the 0.15-2.0 MeV energy range, dependent on 9 - the angle between a line 
connecting the centers of the two detectors and the projection of the direction of the GRB 
from the optical axis in the detector plane - through the Klein-Nishina cross section: 



7 We note that using fa = 1 instead of the relative apparent counts sensitivity, as if only the physical 
apparent cross sections Bij were important, does not change our conclusions. 

8 In fact, the detectors are 8 cm long, comparable to their width and separation. We neglect this width 
throughout the analysis. 



4.2. 



Effect of the Klein-Nishina Cross Section 




(14) 
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k 



1 + ^(1-cos^kn) 



(15) 
(16) 



This will produce directionally preferential scattering from a non-polarized beam of 
gamma-rays in the detectors. If the direction of the GRB is k = sin <px + cos <fiz (where <fi 
is 18deg) the line which joins two detectors as ~r* = cosutx + sin uty, then ^(t) = out and 
$kn(£) — arccos(sin0 cos cut). 

In Fig. 6a, we plot the KN cross section as a function of 7, for 0=18 deg, and for four 
photon energies ranging from 0.15-2.0 MeV, normalized to its maximum value. The cross 
section changes by 20% for photons of energy 0.15 MeV, and by 35% for photons of energy 
2.0 MeV as a detector pair rotates in the detector plane. 

A non-polarized gamma-ray beam will produce an angle- dependent double-event 
count-rate in the RHESSI detectors proportional to the sum of the cross section in two 
detectors separated by 180 deg: 



We show 7kn(t) in Fig- 6b. The KN cross section results in an angle-dependent variation 
in total intensity of 7-9% across the 0.15-2.0 MeV band. For comparison, we include in 
this figure I p (9), for p = 0.045. To within the our observational uncertainties, the KN 
cross-section can always be represented by I p (0) 9 . It is notable that the KN cross-section 
produces an effect on the observed signal which is nearly identical to that expected from an 
intrinsically polarized beam. 

The magnitude of "apparent polarization" (that is, the value of p which parameterizes 
Ikn(0) using I p (8); see Fig. 6c) varies as a function of photon energy; it is p — 5% (for 
photon energy 0.15 MeV), increasing to 5.1% (near 0.30 MeV) and decreasing to 3.7% a 2 
MeV. 

This effect is easy to correct for during data analysis, if one knows the direction of the 
GRB relative to the detector pairs. 

We calculated the fractional magnitude and direction of the KN effect on the 
double-count scattering events, by averaging the magnitude of the effect at the energy of 



9 While the Ikn(0) is not identically parameterizable by I p {9) we find that parameters to I p {9) can typically 
be found which duplicates Jkn(0) to better than 1% for p < 0.1. 



-18- 



each detected count during the GRB (neglecting the 10% of counts due to background) 
and finding the average and standard deviation, which were (pkn) =0.049±0.003. We can 
represent the effect of the KN cross section on the intensity pattern produced by the GRB 
by: 

1 + 0.049 cos(2(# - 9 GRB )) 
/KN " L049 

where #grb is the direction of the GRB projected onto the detector plane (#grb = 45 deg E 
of N). 

In the present data, the magnitude of pkn will be decreased by a factor 
-^scattered /-^observed (< 0.046) (90% confidence). Thus, in the present data, (pkn) < 0.0022, 
which, as we show in §4.3, is well below the detection limit. 

Moreover, the magnitude of this signal is overestimated, as we have assumed infinitely 
thin detectors. The effect of finite sized detectors is to diminish the magnitude of this effect, 
as can be immediately seen, for example, in the limit of infinitely long detectors. Moreover, 
the need for the correction in the present work is absent; the maximum magnitude of 
the effect is 0.049 x830±150=41±7 counts - smaller than the Poisson noise of the total 
double-count events (which is V8230 = 91 counts). Thus, this effect will not be important 
in the present analysis. 



4.3. Results 

In Fig. 5f, we show R{9) for the observation of GRB 021206. First, R(9) is consistent 
with being a constant (reduced chi-squared value x^=1.35 for 11 dof). We therefore find no 
evidence of a signal which might be interpreted as due to polarization in the RHESSI data 
for GRB 021206. 

To obtain a limit on a signal which might be interpreted as due to polarization, we find 
a best-fit for I p (9) (^=1.36, 9 dof) with p < 0.041, or <4.1% (90% confidence), finding a 
(likely, local) minimum at 9 P = 74±22 deg. A similar limit on pkn is obtained holding fixed 
9grb = 45 deg, (pkn < 0.04), which is consistent with our calculated upper-limit for this 
effect. 
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5. Duplication of CB03 Results 

We attempt to duplicate the results of CB03. In particular, we are interested in 
duplicating their observed double-event lightcurve. While the methods of producing the 
double-event lightcurve of CB03 have not been described as of this writing, there are 
only a few variables which were not specified by CB03: (1) AT; (2) whether or not 
multiple-scattering events ("bunches") are retained; and (3) the zero point for 9. 

We show the raw (double-event) counts and residual counts as a function of 9 
(D(9) hereafter) obtained by CB03 in Fig. 7a. The residual counts lightcurve is the 
difference between the raw counts and their MC estimation for the double-event count rate 
(M(9)=D(9)—D nu n(9)), assuming zero intrinsic polarization, and using the GEANT mass 
model for radiative transfer through RHESSI 10 . We note that the average values of D(9) 
and D nul \(9) found by CB03 coincide within <0.1%. 

We obtain a good duplication of the CB03 double-event count curves, using 5T=8 b/is 
(as compared with AT=5 bfis in our previous analysis); retaining all multi-count events 
("bunches", including triples and higher); and setting 9 = to coincide with 12.6 deg E 
of N in the Ecliptic (Fig. 7b, top panel). Here, we find A^ 2 =15540 counts, which is close 
to the 7V 2 =14916 found by CB03. We rescaled the D nuU (9) calculated by CB03 so that 
(-Dnuii(^)) = (D(9)). The residual counts (D(9)—D QU u(9); Fig. 7b, bottom panel) is similar 
to that obtained by CB03 (Fig. 7a, bottom panel). 

Using these same selections, we apply our R{9) analysis. If the modulation seen in 
Fig. 7b (bottom panel) is due to polarization, we expect to see the a sinusoidal modulation 
in R(9), too (see Eqns. 9 & 12). Instead (Fig. 7c), we see that R(9) is scattered, significantly 
different from a constant (x^=2. 49/11 dof; prob.=0.4%), but with no systematic trend 
approximating a sinusoidal modulation. The scatter in R(9) cannot be due to a constant 
polarization signal as implied by CB03 (as we presently show, the scatter is due to the 
"bunches"). We conclude from this that the signal observed by CB03 is not due to 
polarization in GRB 021206. 

We then dropped from the analysis all the "bunches" of counts (multiplicities of 
3 and greater). We do this, because we see no reason that polarized photons should 
be preferentially detected in events which involve >2 scatters. Moreover, as we have 
shown above, the bunches do not follow a Poisson distribution, and do not follow a 



10 It was this MC simulation which gave the estimate of II = 0.19±0.04 - that of all photons which are 
detected as double-count events. 19±4% have not been scattered previously in spacecraft passive material, 
and thus have the polarization intact during measurement. 
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simplistic distribution expected from scattering; we therefore interpret these "bunches" as 
a background signal. There were a total of 1272 such bunches, producing a total of 4592 
double-count counts; after removal, there were 10948 double-count events remaining. We 
performed the analysis again: D{9) and D nn u(9) in Fig. 7d, top panel; with residual counts 
D{9) — D nu n(#)in the bottom panel. While the modulation of the residual counts is reduced, 
there remains a significant deviation (in the range 0=100-180 deg). 

However, the corresponding R{9) (Fig. 7e) is consistent with a constant value (x^=0.66, 
11 dof, prob.=77%). Applying the identical analysis as above, we again obtain a limit of 
p < 0.041 (90% confidence). This means that polarization is not detected (and that the 
scatter in R(8) in Fig. 7c is due to the presence of the "bunches" of counts). 

When we duplicate the analysis of Sec. 2.5 on the data in Fig. 7b and c, we obtain a 
fraction of counts which may be due to double-count scattering events of /=9±2% (x^=6.0, 
6 dof; Prob.=3xl0~ 6 ) and 8±3% (x£=13.9, 6 dof; Prob.=0.03), respectively. The fact that 
data which includes the "bunches" is not well-modeled by our scattering model is consistent 
with the conclusion that these events are not related to scattering. 

To what can the modulation reported in CB03 be due? Since the lack of modulation 
in R(9) implies that D{9) oc S(9), yet D{9) — D nn \\{9) is not constant in 9, the modulation 
must be in D nu n(9), and not in D(9). We attribute, therefore, the modulation observed by 
CB03 to neglected systematic uncertainty in D an u(9), which we explore in § 7. 



6. The Intrinsic Polarization of GRB 021206 

Our parameter p is not the intrinsic fractional polarization II of GRB 021206. The 
relationship of our parameter p to II is: 

p= » n sh < 18) 

where II is the intrinsic fractional polarization of GRB 021206; fi is the fraction of 
i?#£ , 5'S7-detected double-count photons which have not been previously scattered within 
the detector; S is the number of double-count events due to scattering, and B is the number 
of background double-count events (that is, not due to scattering). 

If our parameter p is consistent with zero, then the intrinsic polarization 11 is also 
consistent with zero; thus, regardless of the magnitude of (and uncertainties in) fi, S and 
B, our upper-limit of p < 0.041 implies a non-detection of polarization in the RHESSI data 
of GRB 021206. This is in conflict with the claim of detection described by CB03. 
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We now describe the origin of the discrepancy between these two results. 

Based on our limit p < 0.041, S = 830±150 counts, and S + B = 8230 counts, we find 
the 90% confidence upper-limit for intrinsic polarization in GRB 021206 is II < 0.41//1. 

CB03 finds fi = 0.19±0.04 from their Monte Carlo mass-model simulations, which we 
adopt. Relying upon this value, II < 2.14 (fractional polarization, or <214%) - that is, an 
upper-limit which permits up to 100% intrinsic polarization for GRB 021206, implying that 
the observation is not sensitive to intrinsic polarization in GRB 021206. 

The lack of sensitivity is due to the lack of double-count events from scattering; the 
number we found (830±150 counts) is smaller than and in conflict with the number found 
by CB03 (9840±96). 

If we had found the fraction / of total double-count events from scattering to be that 
found by CB03 (/ =66±1%; see Table 2), we would would instead conclude that II < 0.46 
(90% confidence), compared with the II = 0.80±0.20 (CB03). However, our analysis 
(§2.4) finds that /=11±3%, which produces instead a limit on the intrinsic polarization of 
GRB 021206 well above 100%. 

Our conclusion is that the RHESSI dataset is not capable of detecting polarization 
in GRB 021206, due to the low signal-to- noise of the dataset (that is, a relatively small 
number of double-count scattering events compared to that estimated by CB03). 

7. Systematic Uncertainty in CB03 

In our analysis, the uncertainty in polarization detection is due only to Poisson 
noise (in R(0)); in the analysis of CB03, the uncertainty in polarization detection 
(M(9) = D(6) — D nun (9)) is the sum of Poisson noise and the systematic uncertainty in 
Amii(#)- CB03 neglect the systematic uncertainty in D nu \\(9). Here, we show that this 
systematic uncertainty is not negligible. 

The observable modulation lightcurve is: 

M(9)=D(9)-D nuti (9). 

The systematic uncertainty in M(6) due to radiative transfer through the spacecraft 
mass model appears in D QU \\(9). In a simple formulation, one can represent this as: 
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D aua (6) = (D nun (9)) 

W 

The function and its uncertainty was not given by CB03. To be able to neglect the 
systematic uncertainty in fi, its magnitude would need to be $ half the Poisson uncertainty 
- <20 counts, or 1.6% per angular bin, given the (D(0)) = 1250 counts/bin. 

We can estimate the uncertainty in in the analysis of CB03 from the uncertainty 
in its mean value ((//) = pL — 0.19±0.04), which has a 25% uncertainty. If we take the 
uncertainty in D nu u(8) to be the same fractional value as the uncertainty in the average \x 
(25% per angular bin), we find approximately o" systcma tic ~ 0.25 x 1250(counts/bin) ^310 
counts/bin, much larger than the 20 counts limit, and which should appear directly in their 
difference lightcurve (their Fig. 2, or our Fig. 7a) 11 . However, their analysis shows only 
the Poisson uncertainty (op oisson ~40 counts). Moreover, the neglected £7 S y S tematic = 310 
counts/bin is greater than the reported modulation magnitude (~120 counts/bin); thus, we 
see no reason to neglect it. 

We conclude that the analysis for detection of polarization by CB03 incorrectly 
neglected the dominant uncertainty - the systematic uncertainty D nu u(9), which must 
be added in quadrature with the Poisson uncertainty. This uncertainty, from the above 
approximate analysis, is greater than the reported signal interpreted as due to polarization. 



8. Discussion and Conclusions 

In selecting data to analyze, we found 8230 double-count events (including signal and 
irreducible background). This is well below the 14916 double-count events found by CB03 
(Table 2). The discrepancy appears to be due to two selections employed by CB03: (1) an 
unjustifiably wide time window for "simultaneous" events (AT = 8b/xs, whereas we see only 
AT = 5b/xs as being justified) and (2) the inclusion of multi-count events ("bunches"), in 
which >2 counts are detected within AT. 

We described an analysis to detect polarization in RHESSI data. We apply this 
analysis to the data for GRB 021206, and find no evidence for a signal which might be 
interpreted as due to polarization. The magnitude of polarization-like modulation in the 



11 This is a conservative estimate, as we have assumed the fractional systematic uncertainty in the 12 
individual 9 bins is equal to that of the average value of /u; the fractional uncertainty in the individual bins 
should be larger than that of the average value, by a factor of VT2 if they are Gaussian distributed. 
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lightcurve is p <0.041 (90% confidence). This corresponds to an upper-limit on the intrinsic 
polarization of GRB 021206 of II <214%; that is, we find that the analysis is insensitive to 
polarization at any level in GRB 021206. 

The discrepancy between our derived upper-limit and the claimed detection by CB03 
is due to two effects. First, during detection, CB03 neglected the systematic uncertainty in 
their null lightcurve, due to the Monte Carlo simulation of radiative transfer through the 
RHESSI mass model; the magnitude of this uncertainty is consistent with the magnitude of 
the observed modulation. In comparison, this systematic uncertainty does not play a role 
in the present R(Q) analysis, since the mass model is unnecessary. We are instead limited 
by Poisson noise. Thus, we conclude that the modulation interpreted as due to polarization 
by CB03 is instead due to systematic uncertainty in their null lightcurve. 

Second, during correction of the detection of polarization signals into the magnitude of 
intrinsic polarization in GRB 021206 we find 11±3% of all detected double-count events in 
our selection are due to scattering events. This results in a low signal-to-noise ratio, and a 
resulting high upper-limit on II, such that the observation does not constrain the intrinsic 
polarization of GRB 021206. We justified our data selections and demonstrate that the 
signal should be clean of irreducible background. In contrast, CB03 estimated that 66±1% 
of their detected double-count events were due to scattering events, while not describing or 
justifying their selection criteria. Specifically, our duplication of the analysis of CB03 found 
an unjustifiably larger "simultaneous" window (AT=8 b/is, vs. AT=5 b/is we used) and 
included the multi-count events ( "bunches" ) - both of which increase the background. 

In our analysis, the limiting sensitivity is the highest theoretically possible - set by 
photon counting statistics. We therefore believe our analysis is more robust, and we 
conclude that: (1) the RHESSI observation of GRB 021206 is not sensitive to polarization, 
and (2) there is otherwise no evidence for polarization in the data. 

In conclusion, we find that there is no existing constraint on the intrinsic polarization 
in the gamma-ray flux of GRB 021206 - or, for that matter, any gamma-ray burst. It seems 
unlikely that a constraint will emerge from further RHESSI observations. GRB 021206 was 
extremely bright (brighter than any GRB observed during 8 years of full-sky coverage with 
ULYSSES; Atteia et al. 1999) and detector deadtime already played a role in decreasing 
the signal-to- noise ratio at the peak of the burst, meaning an even brighter burst will 
only moderately improve the signal-to-noise ratio; also, GRB 021206 was located within 
5% of the sky closest to the solar limb, as it must to take advantage of the polarization 
detection approach used herein (distinctly different from that used to measure polarization 
in solar flares McConnell et al. 2002), further decreasing the likelihood of a useful constraint 
on GRB polarization with RHESSI. Such a measurement must therefore wait for more 
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sensitive instrumentation. 
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Fig. 1. — Diagram of RHESSI spectrometer detector layout. Each circle represents the 
7.1 cm diameter detector unit. The numbering of the detectors corresponds to that in general 
use, and these are shown as from the solar perspective. The sense of rotation shown is the 
rotation of the spacecraft; the axis of rotation is not that shown, but is instead close to the 
coordinate origin, (0,0). 

Fig. 2. — (a) Histogram of time-stamp separation of consecutive counts in GRB 021206, in 
units of 2~ 20 s (byus). (b) Same as (a), except consecutive counts in the same detector have 
been replaced by a single count at the earlier timestamp. The peak at non-zero value of AT 
indicates that simultaneous counts can have time-stamps which are not equal. We thus adopt 
the convention that counts with timestamps which differ by < 5 b/is are "simultaneous" . 

Fig. 3.— The 0.15-2.0 MeV lightcurve of GRB021206, as observed with RHESSI, beginning 
at 2002 Dec 6 22:49:10 UT. The two vertical lines delineate the 5.0 s period used in the 
present analysis, which were chosen to match the analysed data in CB03. 

Fig. 4. — Relative frequency of double-count events, by detector. Top panel: Data points, 
with error bars, give the number of double-count events associated with each detector 
and their Poissonian uncertainties. The solid histogram represents the best-fit scattering 
+ coincidence model to these data, which has 11±3% scattering. Dotted histograms 
show the contributions of scattering (lower histogram) and coincidence (upper histogram), 
respectively. Bottom panel: Residuals of the data compared to the model, scaled to the 
Poisson uncertainties. All double-events that interact with detector 2 have been excluded 
from this analysis; see text for details. 

Fig. 5. — In these figures, 9 = and 9 n =0 correspond to celestial north, with increasing value 
rotating east of north about the detector origin (clockwise, facing the Sun; cf. Fig 1). The 
position of GRB 021206 is marked in each panel by a vertical broken line, (a) Single event 
countrate S as a function 9, where 9 is the angle between celestial north and the center of the 
detector in which the count is observed, in the range 0-360 deg. The countrate is inconsistent 
with being constant in 9, with the xl shown, and a probability of observing such variability 
of <10~ 100 . (b) Double-count event countrate D as a function of 9; here, we place 1 count 
for each of the two involved in the event, with each being placed in the appropriate 9 bin. 
The strong variability is inconsistent with being constant in 9, with the given xt, an d a 
probability of of observing such variability of <10~ 50 . (c) The ratio S/D is consistent with 
being constant (prob=0.54). (d) The single-count event rate S in our polarization analysis, 
as a function of sky angle 9 n in the range 0-180 deg (see §4). The countrate is inconsistent 
with being constant in 9 n . (e)The double-count event rate D in our polarization analysis, 
as a function of sky angle 9 n , which has values in the range 0-180 deg (see §4). This figure is 
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directly comparable to the double-count event rate of CB03, without having subtracted the 
Monte Carlo "null" lightcurve. As in CB03, we also find significant variability as a function 
of 9 n . (f) The ratio R = D(9 n ) / S(9 n ) vs. sky angle 9 n . Polarized photons would scatter 
through a preferred angle on the sky, producing more double-count events in that direction 
relative to single count events than in the perpendicular direction, producing variation in 
R(9 n ). However, R is observationally consistent with being constant in 9; thus, we find no 
evidence of polarization. 

Fig. 6. — Panel a: The relative Klein-Nishina cross-section for a pair of detectors as a 
function of angle between the initial photon direction and a line between the two detectors, 
at four different photon energies (listed top to bottom). Panel b: The relative cross section 
for a double-scatter event for a pair of detectors due to the Klein-Nishina cross section, as a 
function of angle between the detectors and the direction of the GRB projected in the plane 
of the detectors, assuming the GRB is offset by 18 deg from the axis of the plane occupied 
by the detectors (zero offset would produce zero modulation). At 18 deg, the magnitude of 
modulation is between 7% (for 0.15 MeV photons) and 9.5% (for 2.0 MeV photons). The 
solid lines are for the same photon energies as given in panel a, with larger energies producing 
greater variation. The dotted line is the 1(9) pattern only for a polarized photon beam, with 
p = 0.045 and 9 P = 0. Panel c: The inferred fractional polarization due only to the K-N 
cross-section, as a function of photon energy. 

Fig. 7. — a) Data copied from Fig. 2 of CB03. (Top Panel) Crosses are the double-count event 
lightcurve D(9) as a function of angle on the sky 9. Solid squares are the MC-simulated null- 
lightcurve D nu \\(9). Note the data is repeated at 180 deg. (Bottom Panel) Residual counts, 
D(9) — D nn n(9), exhibiting the modulation in 9 interpreted by CB03 as due to polarization 
in GRB 021206. The solid-line is their best-fit of a polarization signal (also in panels b 
and d). b) (Top Panel) Our duplication of the results of CB03, finding data selections 
which reproduce their D(9) (see text), with symbols having the same meaning as in panel 
a. The solid points, however, are the D nn \\(9) values calculated by CB03, renormalized (see 
text). (Bottom Panel) The residual counts curve duplicates well that of CB03. c) The 
ratio R(9) = D(9)/ S(9), using our duplicate data selection. The sinusoidal modulation in 
D(9) —D nn n(9) is not present in R(9); instead, the ratio shows significant random scatter, d) 
The same analysis as in panel b, but with the multi-count events removed from the dataset. 
The modulation in D(9)-D nu n(9)is reduced, but still present, e) The ratio R(9) using the 
same data set as panel d. The scatter in panel (c) has disappeared (due to the multi-count 
events removed), and no modulation is present. This implies that D(9) oc S(9) and that the 
modulation in D(9) — Amii(#) seen by CB03 (panel a), and duplicated here in panels b and 
c, is not due to polarization- related modulation in D(9), but in D nvLll (9) (the null lightcurve). 
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Table 1. Detector Physical Parameters 



Detector # 


X 


Y 




h 




(cm) 


(cm) 


(relative solid angle) 


(counts) 


1 


-7.4 


-3.0 


3.73 


9870 


2 


7.4 


-3.0 


n/a 


n/a 


3 


15.4 


5.0 


1.90 


8191 


4 


-15.4 


5.0 


2.68 


8426 


5 


8.2 


14.0 


2.52 


8114 


6 


-8.2 


14.0 


2.94 


7379 


7 


0.0 


5.8 


4.23 


8377 


8 


-5.6 


-14.9 


2.24 


8903 


9 


5.6 


-14.9 


1.80 


7216 


Table 2. Double-Count Events and Irreducible Background 




Events 




Present work 


CB03 


Total Double-Count Events Found 


8230 


14916 


Coincidences 






6640±80 


4488±72 


Other Background" 




760±110 


588±25 


Double-Count Scattering Events' 


•> 


830±150 


9840±96 



Note. — a In the present work this is the background due to the "unknown" mechanism. In 
CB03, this is the value labeled "background scatters". b Using a second, independent method 
which does not require detailed examination of the different background contributions, we 
find 910±250 two-detector scattering events, which lends us confidence in this value. See 
§2.5. 



